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Abstract 

We have searched for the violation of baryon number B and lepton number L in the (B — L)- 
conserving modes r~ — > A°ir~ and r~ — ► A°K~ as well as the (B — L)-violating modes r~ — > 
A 7r~ and r~ -> A if" using 237 ftr 1 of data collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage ring. We do not observe any signal and determine preliminary 
upper limits on the branching fractions B{r~~ — > A°ir~) < 5.9 x 10~ 8 , £>(r~ — > A°7r~) < 5.8 x 10 -8 , 
B(r~ -» A°iTA < 7.2 x 10~ 8 , and B(r _ -> A°K~) < 15 x 10" 8 at 90% confidence level. 
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1 INTRODUCTION 



One of the important unresolved questions of our time is the presence of a large baryon asymmetry 
in today's universe. According to A. Sakharov three conditions must be satisfied in order 
for a baryon asymmetry to arise from an initial state with zero baryon number: baryon number 
violation, C and CP symmetry violation, and a departure from thermal equilibrium. No baryon 
number violating processes have yet been observed |2j. Though we know that the baryon number 
was violated in the early universe we do not know how it came about. Conservation of angular 
momentum requires that the spin 1/2 of a nucleon that is decaying to a lepton be transferred to 
the lepton: AB = ±AL. Therefore there are two types of baryon instabilities \A(B — L)\ = 0,2. 
In the Standard Model (SM), and in most of its extensions, it is required that A(B — L) = 0. The 
second possibility of \A(B—L)\ = 2 allows transitions with AB = —AL, or |AB| = 2 and \AL\ = 0, 
or \AL\ = 2 and \AB\ = 0. It follows that the conservation or violation of (B — L) determines the 
mechanism of baryon instability. 

It has been shown that, in baryogenesis, nonperturbative Standard Model effects at the elec- 
troweak energy scale will erase any baryon excess generated by {B — L)-conserving processes at 
the earliest moments of the universe (T >> 1 TeV) [3j. In addition, generating a baryon excess 
through electroweak effects alone does not seem to be adequate to account for the observed baryon 
asymmetry [I]. A component with A{B — L) = 2 might be necessary to explain baryogenesis. 

Most existing searches for (B—L) violation have been restricted to experiments with nucleons [2]. 
In this analysis we search for the decays r — > An and r —* AK , in the (B — L)-conserving modes 
r~ — > A°tt~ (K~) as well as the (B— L)-violating modes r~ — > A°ir~ (K~). Charge conjugate modes 
are always included if not mentioned otherwise. A similar analysis of the modes r — > Air published 
recently by the Belle Collaboration 5^ finds the upper limits £>(r~ — > yl°7r~) < 14 x 10 -8 and 
B{t~ -> A°tt-) < 7.2 x 1(T 8 at 90% confidence level (C.L.). 

Experimental limits on the proton lifetime imply that the expected branching fraction for r — > 
(p + anything) is not observable in the Standard Model: B(t — > p + X) < 10 -40 The A baryon 
couples weakly to the proton. We would then expect similar but approximately 10 8 times weaker 
[0] constraints from the proton lifetime for r — > Air(K). A recent theoretical paper [2j studied 
dimension-6 operators and concludes that baryon number violation in decays involving higher 
generations, assuming proton stability, will not be observable. However such a model may not be 
adequate to describe the apparent baryon asymmetry in the first place. Models with dimension-9 
operators and yet unknown mechanisms that generate baryon number violation or enhance the 
coupling to higher generations may be able to accomplish this (Sj. 

With the advent of the B factories, that also produce large quantities of r leptons, we are now 
able to experimentally study such decays with greatly improved precision. 

2 THE BABAR DETECTOR AND DATASET 

This measurement was performed using data collected by the BABAR detector at the PEP-II storage 
ring. Charged particles are detected and their momenta measured by a combination of a silicon 
vertex tracker (SVT), consisting of 5 layers of double-sided detectors, and a 40-layer central drift 
chamber (DCH), both operating in a 1.5-T axial magnetic field. Charged particle identification is 
provided by the energy loss in the tracking devices and by the measured Cherenkov angle from an 
internally reflecting ring-imaging Cherenkov detector (DIRC) covering the central region. Photons 
and electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The EMC is sur- 



8 



rounded by an instrumented flux return (IFR). Electrons are identified using measurements from 
the DCH, EMC, and DIRC. The average identification efficiency is approximately 97%, whereas the 
pion (kaon) misidentification rate is less than 2% (1%). Kaons are identified using the SVT, DCH, 
and DIRC. The average identification efficiency for the tight kaon selection is approximately 80%, 
whereas the pion misidentification rate is less than 1%. The average identification efficiency for the 
loose kaon selection is approximately 90%, whereas the pion misidentification rate is less than 4%. 
Protons are identified with a likelihood based algorithm using measurements from all described 
detector components. The proton identification efficiency ranges from approximately 90% to 96% 
depending on polar angle and momentum, whereas the average pion (kaon) misidentification rate 
is 5% (12%). Details of the detector are described elsewhere [2]. 

The data sample used corresponds to an integrated luminosity of 237 fb^ 1 collected from e + e~ 
collisions at, or 40 MeV below, the T(AS) resonance. Production and decay of the tau leptons are 
simulated with the kk2f jlOllllj and tauola |12| ll.'Sj Monte Carlo (MC) event generators, according 
to two-body phase space, and taking spin correlations into account for the signal mode. B meson 
decays are simulated with the EvtGen generator ^1], and qq events, where q = u,d,s, or c quark, 
with the JETSET ^)] generator. The detector is fully modelled using the GEANT4 simulation 
package [Tfij . 

3 ANALYSIS METHOD 

We reconstruct candidate events e + e _ —* t + t~ with one r decaying to Air{K) and A — > pn. The 
other tau in each event is required to be a one-prong decay. Decays that conserve (B — L) are 
recognized by opposite sign charge of the pion or kaon from the r decay and the pion from the A 
decay. In decays where (B — L) is violated the two charges have the same sign. 

Each event must have exactly four well reconstructed tracks in the fiducial volume of the DCH 
with a total charge of zero. We divide the events into two hemispheres defined by the thrust axis 
of the event. The thrust axis is calculated using tracks in the drift chamber and calorimeter energy 
depositions without an associated track. We require that the three signal tracks are contained in 
one hemisphere and that there is exactly one remaining track in the other hemisphere, which we 
will refer to as the tagging hemisphere. 

One of the signal tracks must be identified as a proton and, when combined with an oppositely 
charged signal track, must give a p-K~~ invariant mass within 5 MeV/c 2 of the nominal A mass 2 . 
The set of signal tracks are subjected to a topological fit to the decay tree r —* Air(K), which must 
converge and return a x 2 probability greater than 2.5%. 

We require that the center-of-mass (CM) momentum of the A° is greater than the lower kine- 
matic limit of 1.8 GeV/c for r~ — > A°ir~ decays. A requirement on the A flight distance L^o > 1 cm 
and the signed flight length significance L^/a^o > removes t + t~ (88%) and qq (22%) events 
that do not contain true A particles. The remaining backgrounds are mostly from qq events and 
to a lesser degree t + t~ events that contain K® decays and photon conversions 7 — ► e + e~ . None of 
approximately 800 million MC BB events survive the selection criteria. 

Figure n shows a comparison of the MC simulation with our data. Note that the A momentum 
spectrum shown in Figure Q ( a ; b) is not very well described by our MC simulation. This is most 
likely due to imperfections of the qq MC event generator. For this reason the final background will 
be determined from the data. All other variables that were studied show better agreement between 
data and MC. 

We require that the pion track from the A decay as well as the tagging track from the other 
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Figure 1: The A candidate momentum in the e + e~ rest frame for (a) r — ► An and (b) r — > yli^T 
and the yi° invariant mass spectrum for (c) r — > yl7r and (d) r — ► yli^T. The solid stacked histograms 
are from top to bottom: uds backgrounds (blue), and t + t~ (green). The cc component is too small 
to be seen in these figures. The signal Monte Carlo distributions are shown with the dashed red 
histogram. The points correspond to events in the data sidebands. 



r lepton do not pass tight kaon identification requirements. In the mode r — ► Air we require 
that the ir is not identified as a kaon. In the mode r — > AK we require that the kaon track be 
identified with loose kaon identification requirements. To suppress candidates that include tracks 
from photon conversions, we require that the pion or kaon from the r decay and the pion from the 
A decay must not be identified as an electron. The pion or kaon from the r decay must not be 
identified as a proton. 
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We study events in the two dimensional plane m^tx) versus AE An ( K ), where m^ w ^x) is the 
invariant mass of the A and the pion (or kaon) candidate, and AE^i K \ = E An ^ K ) — yfs/2 is the 
reconstructed energy E An r K ) of the signal tracks minus the expected r energy, which is half the 
known e + e~ center-of-mass energy yfs. A rectangular region that includes the signal region was 
blinded during the development of this analysis. Signal candidates are counted in an elliptical 
signal region with a half width of 10 MeV in m A7T ^x) and 90 MeV in AEfafK) centered around 
the nominal r mass and AEj^r K ) = 0. In the case of r — > AK the width in m^rx) is reduced 
to 7 MeV because of the better resolution in this mode. The elliptical signal region is slightly 
tilted to reflect the small correlation between the two variables. The tilt is ~ 3°, which can also be 
expressed as a correlation coefficient between the two variables: p = 0.42 for r — ► Air and p = 0.56 
for r — > AK. The definition of the signal region as well as the other selection requirements applied 
in this analysis have been optimized using MC simulation, to obtain the lowest average upper limit 
for the signal modes under the assumption that no signal will be observed. 

We estimate the number of background events in the signal region with a 2D unbinned maximum 
likelihood fit of the m^rx) an d ^Ej^ik) distributions outside the blinded region. We try a 
number of functional forms that describe both the data and MC distributions. The default fit uses 
a simple parametrization that describes the data well and results in a background estimate that 
is in the center of the possible range of values. A first-order polynomial is fitted to the m^^x) 
distribution and a Gaussian function to the AE^rx) distribution. The blinded region is excluded 
from the fit and the probability density function is set to zero within the blinded region. The 
parametrizations obtained are shown in Figure |2j The elliptical signal regions and the blinded 
region are also indicated in Figure |31 Due to the uncertainties of the background parametrization 
and the possibility of correlations among the fit variables, we take a conservative 100% error on 
the number of estimated background events in the signal region. 

4 SELECTION EFFICIENCY 

The signal efficiencies have been obtained from Monte Carlo simulations. Systematic uncertainties 
have been studied using independent control samples of real data; a summary is presented in 
Table ^ The largest contributions are from uncertainties related to the tracking efficiency, and 
A reconstruction. The latter has been estimated by comparing lifetime distributions of long lived 
particles in data and Monte Carlo. The uncertainty on the branching fraction B(A° — > pir~) has 
been taken from the Review of Particle Physics 0. Contributions to the systematic uncertainty 
are added in quadrature to give a total systematic uncertainty of 6.9% in the mode r — > Air and 
7.0% for r -» AK. 

5 RESULTS 

The data distributions in the AEj^rx) versus m^iK) plane after all selection requirements are 
shown in Figure No signal candidate events are observed in the r — * Air mode. We observe 
one candidate event in the (B — L)-violating mode r~ — ► A°K~. We determine upper limits on 
branching fractions at 90% C.L. using the method described in Ref. |17|. This method considers 
uncertainties both on the signal efficiency as well as the number of expected background events in 
the signal region. The number of expected background events and number of observed events in 
the signal region, the signal efficiency, and the upper limit that has been determined are shown 
separately for the (B — L)-violating and {B — L)-conserving cases in Table El The upper limit on 
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Figure 2: Projections of the background parametrization as derived from a 2D unbinned maximum 
likelihood fit. The top row shows the AE projection, and the bottom row the rriAh projection. The 
r — ► Air mode is shown in the left column and the r — ► AK mode in the right column. The fitted 
probability density function (PDF) is indicated by a line. The PDF was required to be zero in the 
blinded region, which causes the apparent drop around the signal regions in these projections. The 
points with error bars correspond to r — ► Air(K) candidates in data, outside the blinded region. 



the branching fraction is given by 

where t is the 90% C.L. upper limit for the signal yield, a TT = 0.89 nb is the assumed cross section 
for production of r pairs, C = 237 fb _1 is the total luminosity of our dataset, B(A — > pn) = 0.639 
is the A branching fraction taken from the RPP J2], and e is the signal efficiency. 
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Table 1: Summary of systematic uncertainties on the signal efficiency and the luminosity and cross 
section. 



suLii ce 


UIlCcI IdllJXy \ /0 ) 


A reconstruction 


5.0 


tracking efficiency 


4.0 


proton identification 


1.0 


kaon identification (r — > AK only) 


1.0 


B{A -»• pir) 


0.8 


luminosity and cross section 


2.3 


total r — > yl7r 


6.9 


total r — »■ ylif 


7.0 



Table 2: The number of expected background events in the signal region, signal efficiency, number 
of observed events, 90% C.L. upper limit for the signal yield (£), and the upper limit branching 
fraction for each mode. 



mode (B — L) expected efficiency observed £ upper limit on B 









background 


% 


events 




@ 90% C.L. 




A°ir- 


conserving 


0.42 ± 0.42 


12.28 





1.97 


5.9 x 10^ 8 


T~ 


-^A°TT- 


violating 


0.56 ±0.56 


12.21 





1.90 


5.8 x 10" 8 


T~ 


-»• A°K- 


conserving 


0.26 ±0.26 


10.63 





2.08 


7.2 x 10" 8 




—> A°K~ 


violating 


0.12 ±0.12 


9.47 


1 


3.78 


15 x 10~ 8 



6 SUMMARY 

A search for the (B — L)-conserving modes t~ — ► yl°7r~ and t~ — ► A°K~ as well as the (B — L)- 
violating modes — > 7l 7r~ and r~ — ► A°K~ has been performed using 237 fb _1 of e + e~ data. 
No signal is observed and we obtain preliminary upper limits on the branching fractions at 90% 
C.L. of B{t~ -»• A°7r-) < 5.9 x 10~ 8 , B(t~ -► /l ^-) < 5.8 x 10~ 8 , jB(t" -► 4°^-) < 7.2 x 10~ 8 , 
and B(t~ — ► yl ivT _ ) < 15 x 10~ 8 . This analysis is the first measurement of the mode r — ► AfsT, 
and it improves over earlier measurements of the mode r — ► yLr. 
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Figure 3: AE^ n r K \ versus ttij^ik) data distributions for the (B — L)-conserving modes (left) and 
the (B — L)-violating modes (right). The top row shows the mode r — ► Air; the mode r — ► AK 
is shown in the bottom row. The expected signal distribution (taken from Monte Carlo) is shown 
with red squares; data events are shown as dots. The large rectangles in each plot are from left to 
right: left sideband, blinded region, and right sideband. The elliptical signal region is also shown. 
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